Simulation of fault-tolerant quantum circuits on quantum computational tensor 

network 



Tomoyuki Morimae 1,2 and Keisuke Fujii 3 

1 LAMA, Universite Paris-Est Marne-la-Vallee, 77454 Marne-la-Vallee Cedex 2, France 

2 IRCS, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8550, Japan 

3 Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan 

(Dated: February 16, 2012) 

In the framework of quantum computational tensor network [D. Gross and J. Eisert, Phys. Rev. 
Lett. 98, 220503 (2007)], which is a general framework of measurement-based quantum computa- 
tion, the resource many-body state is represented in a tensor-network form (or a matrix-product 
1 form), and universal quantum computation is performed in a virtual linear space, which is called 

a correlation space, where tensors live. Since any unitary operation, state preparation, and the 
projection measurement in the computational basis can be simulated in a correlation space, it is 
£Sj ' natural to expect that fault-tolerant quantum circuits can also be simulated in a correlation space. 

However, we point out that not all physical errors on physical qudits appear as linear completely- 
positive trace-preserving errors in a correlation space. Since the theories of fault-tolerant quantum 
circuits known so far assume such noises, this means that the simulation of fault-tolerant quantum 
circuits in a correlation space is not so straightforward for general resource states. 
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I. INTRODUCTION 



Quantum many-body states, which have long been central research objects in condensed matter physics, statistical 
physics, and quantum chemistry, are now attracting the renewed interest in quantum information science as fun- 
damental resources for quantum information processing. One of the most celebrated examples is one-way quantum 
computation [H-Q ■ Once the highly-entangled many-body state which is called the cluster state is prepared, universal 
quantum computation is possible with adaptive local measurements on each qubit. Recently, the concept of quan- 
tum computational tensor network (QCTN) 04^]) which is the general framework of measurement- based quantum 
computation on quantum many-body states, was proposed. This novel framework has enabled us to understand how 
general measurement-based quantum computation is performed on many other resource states beyond the cluster 
state. The most innovative feature of QCTN is that the resource state is represented in a tensor-network form (or a 
matrix-product form) [7j-(9j , and universal quantum computation is performed in the virtual linear space where tensors 
live. For example, let us consider the one-dimensional open-boundary chain of N qudits in the matrix-product form 

mL,R))? = ■= = = E - E (L\A[k N }...A[k 1 ]\R)\k N ,..,h), (1) 
VJn\\ l /> \ R I) fe 1= o 



k N =Q 

where 

f N (\L),\R)) = {L\(A N \R){R\)\L) 

is the normalization factor, 



i=0 

is a map, {|0), \d — 1)} is a certain basis in the d-dimensional Hilbert space (2 < d < oo), \L) and \R) are D- 
dimensional complex vectors, and {j4[0], A[d — 1]} are D x D complex matrices. Let us also define the projection 
measurement Me,4> on a single physical qudit by 

M e ,4 = {M,\M,\2),...,\d-l)}, (2) 

where 

K ) = cos||Q) + e^sin||l), 
|j8 M > = sin^|0)-e^cos^|l), 
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< 9 < it, and < <fi < 2ir. If we do the measurement M.e,<t> on the first physical qudit of \^(L, R))i and if the first 
physical qudit is projected onto, for example, \otg </,) as a result of this measurement, the state \^(L, R))i becomes 
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d-l d-l 



— , ir ., —r^ J2 - £ (L\A[k N }...A[k 2 }A[ae,t\\R}\k N ,...,k 2 }®\a e ,t) 

VfN-i(\L),A[ae,4,]}R))j^ fcjy=0 



VfN-i(\L),A[a e J\R)) t± ^ \\A[ae,i 



-\R)\k N ,...,k 2 ) ® \a e><t> ), 



where 



Then, we say "the operation 



A[ag^} = cos -A[0] + e~ i0 sin -A[l] 



is implemented in the correlation space". In particular, if A[0], A[l], 9, and <f> are appropriately chosen in such a way 
that A[o!0,0] is proportional to a unitary, we can "simulate" the unitary evolution 



\\AM\\ 



of the vector \R) in the virtual linear space where A's, \R), and \L) live. This virtual linear space is called the 
correlation space The core of QCTN is this "virtual quantum computation" in the correlation space. If 

the correlation space has a sufficient structure and if A's, IL), and \R) are appropriately chosen, we can "simulate" 
universal quantum circuit in the correlation space 0-0, Il0l4l3 | . 

For the realization of a scalable quantum computer, a theory of fault-tolerant (FT) quantum computation [R 
necessary. In fact, several researches have been performed on FT quantum computation in the one-way model 
I22I ]. However, there has been no result about a theory of FT quantum computation on general QCTN [2c 
particular, there is severe lack of knowledge about FT quantum computation on resource states with d > 3. 
necessary to consider resource states with d > 3 if we want to enjoy the cooling preparation of a resource state and 
the energy-gap protection of measurement-based quantum computation with a physically natural Hamiltonian, since 
no genuinely entangled qubit state can be the unique ground state of a two-body frustration- free Hamiltonian (28| . 

One straightforward way of implementing FT quantum computation on QCTN is to encode physical qudits with a 
quantum error correcting code: 




d-l 



|*) = ) £ ... £ (L\A[k N ]...A[k 1 ]\R)\k N ,...,k 1 ) 
Vjn(l,R) kl=0 



k N =0 



where \ki) (i — 1, ... JV) is the encoded version of (such as |0) = |000) and |1) = | 111) , etc.) In fact, this strategy 
was taken in Refs. |2ll [22| for the one-way model (d = 2), and it was shown there that a FT construction of the 
encoded cluster state is possible. For d > 3, however, such a strategy is difficult, since theories of quantum error 
correcting codes and FT preparations of the encoded resource state |^) are less developed for d > 3. Furthermore, if 
we encode physical qudits with a quantum error correcting code, the parent Hamiltonian should no longer be two-body 
interacting one. 

The other way of implementing FT quantum computation on QCTN is to simulate FT quantum circuits in the cor- 
relation space. Since any unitary operation, state preparation, and the projective measurement in the computational 
basis can be simulated in a correlation space (for a more precise discussion about the possibility of the measurement, 
see Ref. [13]), it is natural to expect that FT quantum circuits can also be simulated in a correlation space. An 
advantage of this strategy is that theories of FT quantum circuits for qubit systems are well developed |13l - fl7| . In 
fact, this strategy was taken in Refs. [H, [l9| for the one-way model (d = 2). They introduced a method (which 
we call "the ensemble method" since the ensemble of all measurement results arc considered) of simulating quantum 
circuits in the correlation space of the cluster state, and showed that all physical errors on physical qubits can be 
linear completely-positive trace-preserving (CPTP) maps in the correlation space of the cluster state. This means 
that FT quantum circuits can be simulated in the correlation space of the cluster state. 
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In this paper, however, we point out that it is not so straightforward to simulate FT quantum circuits in a correlation 
space of a general resource state. In the next section, Sec. [TT1 we review the simulation of FT quantum circuits on the 
one-dimensional cluster state [HI, in terms of the QCTN picture to fix the notation. We see that for the cluster 
state all physical errors can be linear completely-positive trace-preserving (CPTP) maps in the correlation space, and 
therefore the theory of FT quantum circuits can be used in the correlation space. However, this is not the case for 
other general resource states of QCTN. As an example, we consider a similar way of simulating quantum circuits in 
the correlation space of the one-dimensional AKLT state [ID, H3] in Sec. Mil and show that not all physical errors can 
be linear CPTP maps in the correlation space of the AKLT state. Since all theories of FT quantum circuits known 
so far assume such noises [l3T - fr7j , this means that it is not so straightforward to apply these FT theories to quantum 
circuits simulated in the correlation space of general QCTN. In Sec. IIV1 we give some intuitive explanations of the 
reason why the cluster state is so special, and why not all resource states work as the cluster state. In Sec. [V] we 
consider another standard way of simulating quantum circuits in the correlation space, which we call "the trajectory 
method" since a specific trajectory (measurement results) is considered. However, we show a general theorem that 
such an another way does neither work if d > 3. 

In short, we show in this paper that it is not so straightforward to simulate FT quantum circuits in the correlation 
space of a general resource state. Since all errors behave nicely in the correlation space of the cluster state [laUJii l ess 
attention has been paid to the difference between a real physical space and a correlation space of a general resource 
state. Our results here suggest that these two spaces can be different, and because of the difference, simulations of FT 
quantum circuits can be difficult in a correlation space. Of course, we do not show here the impossibility of making a 
QCTN fault-tolerant. In a future, a highly elaborated method might be found which makes all QCTN fault-tolerant. 
We hope that our results will help to study such a challenging subject of a future study. 

Assumptions: Throughout this paper, we make the following assumptions: Since the MPS R))i , Eq. ([1]), is 

a resource state for measurement-based quantum computation, we can assume without loss of generality that Afa^], 
A[/30 j( ji], A [2], A [3], A[d — 1] are unitary up to constants: 



A{ a e,<t>] 




A\fie,4>] 


= cpUp, 


A{2] 


= c 2 U 2 , 


A[3] 


= C3U3, 


A[d-1] 


= Cd-iUd 



(3) 

where c a , cp, c 2 , ■■■ Cd-i are real positive numbers, U a , Up, U 2 , Ud-i are unitary operators, and 

A[(3 g J = sin 6 -A[Q] - cos ^A[l]. 

This means that any operation implemented in the correlation space by the measurement M.e,<j> on a single physical 
qudit of \^(L, R))^ is unitary. Note that this assumption is reasonable, since otherwise \^f(L,R))^ does not seem 
to be useful as a resource for measurement-based quantum computation. In fact, all known resource states so far [l|- 
[29|, [3ll - l33l | , including the cluster state and the AKLT state, satisfy this assumption by appropriately rotating each 
local physical basis. Furthermore, we can take c a , cp, c 2 ,..., Cd-i such that 

d-l 
k=2 



since 



(L\A[k N ]...A[k 1 ]\R)\k N ,...,k 1 ) = =L== VC N Yl (L\^...^\R)\k N ,...,h) 

u ..., kN VM\ L )AR)) kl ,...,k N vc vc 



VfN(\L),\R)) kl ^ kN ^f N {\L),\R)) 
and we can redefine A[ki]/\f(J — > A[ki 



II. SIMULATION ON THE CLUSTER STATE 



Let us first review the results for the cluster state [H, El in terms of the QCTN picture to fix the notation. 
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A. Simulation on the cluster state without error 



Let us first assume that there is no error. The one-dimensional cluster state is the matrix-product state defined by 
d = 2, A[0] = |+)(0|, and A[l] = | — We measure each physical qubit in the basis 

\6.) = ^=(|0> + (-l)»e*|l», 

where s e {0,1}. In the correlation space, X S J(6), where J{6) = He MZ / 2 , is implemented. First, we measure the 
first physical qubit in the {\6q), \0i}} basis. Then we obtain 



2f N (\L),\R)) 



J2 W(X s iJ(6)\R)) 2 ®\0 Sl }(6 Sl \® m ( Sl ), 



where m(0) and m(l) are mutually orthogonal states which record the measurement result, and 

W{\i,)) r = E (i|^Ar]...A[A ;r ]|^)(^|At[fc;]...At[^]|L)|^,...,M<^,-",fc;i- 



k r ,...,kjy k' r ,...,k' N 

If we trace out the measured first physical qubit \0 Sl ){0 Sl |, we obtain 



si=0 

Second, we measure the second physical qubit in the {X Sl \<f>o), X Sl \<f>i)} basis. Then we obtain 
i 1 1 

92f nn rsu E E W(X»J((-iy*4>)X* J(6)\R)) 3 ® X s ^ S2 ){<t> S2 \X^ ® m( Sl ) ® m(s 2 ) 
2 f N {\L),\R)) s ^ Qs ^ 

1 1 1 

= 92 f fin rm^ E E ^(^ S2 ^ S1 J(f?)| J R» 3 ® X- |^ 2 )<^ 2 \X'* ® m( Sl ) ® m( S2 ). 

If we trace out the measured second physical qubit X Sl \(j) S2 )((j) S2 \X Sl , we obtain 

1 1 1 

2*f N QL),\R)) S E W(^ S2 ^ S1 ^)^)| J R))3 ® m( Sl ) ® m(s 2 ). 

Third, we measure the third physical qubit in the {Z Sl X S2 \rjo), Z Sl X S2 \i]i)} basis. Then we obtain 

ill 

E E E mx S3 x si J((-n s2 v)x s2 z s ij^)j(d)\R)) 4 ® z^x^x^ix^ 



2 3 /jv(|i),l^)) 

J« vi / I // S3=0s2= o Sl =o 

®m(si) ® m(s 2 ) ® m(s 3 ) 
iii 

E E E W^' 8 ^ 3 ■/(»?) J(0)J(6»)|i2)) 4 ® Z'-X"\r, aa )(ri 8a \X"Z^ 



2 3 f N (\L)AR)) 

/ I // S3= o S2= o Sl =0 

®m(si) ® m(s 2 ) ® m(ss). 
If we trace out the measured third physical qubit Z S1 X S2 \r] S3 ) (rj S3 \X S2 Z S1 , we obtain 

1 iii 

93 f nr\ \B\\ E E E W(X««Z'V(^) J(<?)|fl)) 4 ® m( Sl ) ® m( S2 ) ® m( S3 ) 

i 1 1 

= 93 f nr\ E E ^(^ S3 ^V(ry) J(0) J(0)|i?)) 4 ® / ® m(s 2 ) ® m( S3 ). 

In this way, we can simulate the desired unitary operation J(r))J(<f))J(6) on \R) up to Pauli byproducts X S3 Z S2 in 
the correlation space. These Pauli byproducts can be corrected later, since they are specified by m(s 2 ) ® m(ss). 



B. Effect of a CPTP error on a physical qudit of a general state 

Before studying the simulation of quantum circuits on the cluster state with error, let us consider the effect of a 
CPTP error on a physical qudit of general resource states, since we will use it later. Let us assume that a CPTP error 

w 

where Y^j=\ E j E j = ^> occurs on the first physical qudit of \^(L, R))i'. 

fN QL) \R)) E E E <i|^[fcjv]--.A[/ Sl ]|J2><i2|^t[*^]...>l+[A^]|i>|A: JV ,...,A ;2 ><A^ ® i^lAsiXAi . 

j k± .... .kjy k'^ .... , k' N 

If we measure the first physical qudit in a certain basis {|m s )}, 

j 5 s k i , . . . , k at k'^ , . . . -,k'jy 

®\m s ){m s \F j \k 1 ){k' 1 \F^\m s ){m s \ 
= f \rY) E E E <i|>l[fciv].-.A[^]^, s |i2><i2|St B A+[^]...At[fe^]| J L>|fc JV fcaXfci, ® Im^Xm.l 

= ^p^yE(E^>l^))>K)^l' 

where 

E J,s = ^2 A [ k K m s\Fj\k)- 
k 

If we trace out \m s ), we obtain 

T^EE-^ 

- 7 jT^ m E E <i|^[fej>r].-.^[fe](S^.-l^><^l^,-) At [^]-" At [^l i >l^'"-' fe 2><^ A;^|. 

• Wl 1 " k 2 ,-,k N k' 2 ,...,k' N j,s 

This means that the map 

\R)(R\^J2 E ^( R \ E l ( 4 ) 

is implemented in the correlation space. Note that 

Y, E l E J,s = E ^[k]A[k']{k\F}\m s )(m s \F 3 \k') 

j,s j,s,k,k' 

= ^A^[k]A[k']{k\k') 

k,k' 

= Y,Ai[k]A[k>](k\U M ,,*ul lt Jk>) 

k,k' 

d-1 

= A^[a e ^]A[a e ,^+ A^[P8^]A[P8^] + ^A^[k]A[k] 

fe=2 

= I, 
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where 

d-l 

u Me ,, = M(o\ + \Po,*){i\ + 53 1^1 

k=2 

is a unitary operator. Therefore, the map Eq. (Q} is CPTP. 

Note that this result does not mean that we can always have CPTP errors in the correlation space: In this section, 
we did not consider any quantum gate. If we implement quantum gates in the correlation space, the situation becomes 
more complicated, and, as we will see later, we sometimes have non-CPTP errors in the correlation space. 



C. Simulation on the cluster state with error 



Now let us consider the case where we implement quantum gates on the cluster state with error. We assume that 
a CPTP error occurs on the first physical qubit of the one-dimensional cluster state. If we measure the first physical 
qubit in the {|#o)j |$i)} basis after such an error, we obtain 

, 1 

f N n L \\ R)) E ( E W ( E ^ \ R )h) ® W«i I ® m ( fi i)- 

si — j 

By tracing out the first physical qubit \9 Sl )(0 Sl \ , we obtain 

£ ^ ? W[EhSl m ^ ® W(Sl) ' 

Second, we measure the second physical qubit in the {X Sl |0 O ), X Sl \<j>i)} basis. Then we obtain 
1 1 1 

2f n L) ]B)) E E ( E W(X S2 J((-l) s ^)E j>Sl |i?)) 3 ) ® X s > \^ 2 )(<j> S2 |X S1 ® m( Sl ) ® m( S2 ). 

jn[\ /, i /; S2=0si=0 j 

By tracing out the second physical qubit X Sl \<j) S2 )((f) a2 \X Sl , 

1 1 1 

of TTn imS E E (E^^ 2 ^^ 1 ) 81 ^^!^)^ ®rn( Sl )®m(s 2 ). 
zW),\n)) S2=0si=0 j 

Third, we measure the third physical qubit in the {Z Sl X S2 \i]q) , Z Sl X S2 \r]i)} basis. Then we obtain 

■'-iVU hi II S3lS2jSl j 
®m(si) S3 to(s 2 ) ® m(s3) 
1 



J2 (^W(X s «X^Z s V(r ? )J((-l) s ^)^ >1 |E})4)®^ 1 X^| % 3><r 7s3 |X^^ 



(g)m(si) ® m(s 2 ) ® m(s3). 
If we trace out the third physical qubit Z Sl X S2 \r) S3 ){r) S3 \X S2 Z Sl 
1 



2 2/Ar (|L),|i?» ss . ^ 

If we further trace out the first record state ra(si) 
1 



53 (E^^ 153 ^ 81 ^ 2 " 7 ^)" 7 ^" 1 ) 151 ^^ 8 !!^)) 4 ) ® m(si) ® m(s 2 ) O m(s 3 ). 



2 2 /iv(|i),|i?)) 



53 (VpP 3 Z S2 J(r?) |i?)) 4 + W(X S3 XZ S2 J(r?) J(-0)%i|i?)) 4 ) <8 m(s 2 ) ® m(s 3 ) 



1 



2 2 M|i),|i?)) S; 



53 (V(X«Z S2 JO?) J{(f>)E jtQ \R)) 4 + W(X S:i Z S2 J(rj) J{<j))XE hl \R))^ ® m(s 2 ) <g> m(s 3 ). 
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In other words, the map 

\ R H R \ E (%o|i2)<i«l4o +^i,ilWl4i*) 
J 

is implemented in the correlation space up to the rotation J(r))J(<fr) and a Pauli byproduct X S3 Z S2 . 
Since 

E (4o%o + (4 1 x)(x^, 1 )) = E4«^> = 

3 ],s 

this is CPTP error. In short, a CPTP error on a physical qubit becomes a CPTP error in the correlation space of 
the cluster state. As is shown in Appendices lAl and [Bl similar result is obtained for the tricluster state [32j . which is 
a variant of the cluster state. 



III. AKLT STATE 



We have seen in the previous section that CPTP errors on physical qubits of the one-dimensional cluster state 
become linear CPTP maps in the correlation space. However, this is not always the case for general resource states. 
In order to see it, let us consider the one-dimensional AKLT state as an example. 



A. Simulation on the AKLT state without error 



First we assume there is no error. The one-dimensional AKLT state is the matrix-product state defined by d = 3, 

m = 

m = -Lxz, 

If we measure the first physical qutrit of the AKLT state in the basis 

( 6 9 6 ^ 

= { cos - 1 0} + i sin - 1 1) , sin - 10) - i cos - 1 1} , |2) ] , 

we obtain 

1 2 

~, ,, n , p n E w (Qi( s i)\ R )h ® Pi{*i) ® m(«i), 
Sf N {\L), \H)) ^ Q 

where pi(si) is the state of the first physical qutrit after the measurement, m(si) is the register state which records 
the first measurement result, and 

Qi(0) = XS Z (6), 
Qi(l) = XZSz(0), 
Qi(2) = Z. 

Here, Sz{9) = e~ lZB l 2 . By tracing out the first measured physical qutrit pi(si), we obtain 

Next we measure the second physical qutrit by choosing measurement basis according to si [2{|. Then, we obtain 

f 2 2 

q2 , Tfn inw E E ^(^(si, s 2 )Qi(si)|i?)) 3 ® p 2 (s 2 , si) ® m(si) ® m(s 2 ), 



8 



where p2(si,S2) is the state of the second physical qutrit after the measurement, m(s 2 ) is the register state which 
records the second measurement result, and 

Q 2 («i,0) - XSz(0), 
Q 2 («i,l) = XZS Z {6), 
Q 2 (si,2) = Z, 



if si = 2 and 

Q 2 (*i,0) = X, 
Q 2 («i,l) = XZ, 
Q 2 (*i,2) - Z, 

for other si. By tracing out the measured second physical qutrit p 2 (s 2 , Si), 

j 2 2 

If we repeat these process, after measuring the rth physical qutrit, we obtain 

^ 2 2 

o rf 7TT\ im^ Z ••• Z W(<9r(si, ■ •' Sr)-Q 2 (si, s 2 )Qi(si)|i?)) r+ i <g> ra(si) <g> ... <g> m(s r ), 
3 ^ 



where 



Qfc(«i,...,s fc _i,0) = XSz(0), 



if si = ... = Sk-i = 2 and 



Qk(si, 


Sfc-1,1) = 


*ZS Z (#), 


Qk(Sl, 


Sk_i,2) = 


z, 


Qk(Sl; 


...,Sfe_i,0) 


= X, 


Qk{Sl; 


— ,Sfe-l, 1) 


= XZ, 


Qk(si, 


...,Sfe_i,2) 


= z, 



for other si, Sfc_i. 
If si = ... = s r = 2, 



For other s\, s r 



where 



Q r (si, s r )...Q 2 (si, s 2 )Qi(si) = Z r . 

Q r ( 8l , 8 r )...Q 2 (8 1 ,8 2 )Ql(8 1 ) = Z^ S Z (0), 



/(si,...,s r ) = 0(5 SiiO 8 5 Si ,i), 

i=l 

r 

5(si,...,s r ) = ®(tf Si ,i e 5 Si;2 ). 

i=l 

Let us add the flag state ij as 

o rf 77T\ i D n Z ^(^(si, Sr) -Q2 (si , « 2 )<3l (si ) | R))r+1 ® m(si) <g> ... <g> m(s r ) <g> ??(/(si, S r ),0(si, S r )) 

3 jN\m, i-K/j „, „ 



9 



where 77(0, 0), 77(0, 1), 77(1, 0), and 77(1, 1), are mutually orthogonal with each other. If we trace out all register states, 
m(si), m(s 2 ), and m(s r ), we obtain 

0r r 7|7\ i p \n E W{Qr(si,...,Sr)-.Q2(si,S 2 )Ql(si)\R))r+l<E)r](f(s 1 ,...,Sr),g(si,...,S r )) (5) 

a r JN{\L), \H)) J— 1 



si,...,s 

1 1 



where 



and 



~ rf ( L E E (KJ ' ^(^^^WI^Wi + MP, 9, r)W(iH#) Wi) ® r?(F, g), (6) 



S r p ^ q = {fa,..., s r ) e {0, 1, 2} xr \ (2, 2) | / ( Sl , 8r ) = p and g(s u s r ) = q} 



h(v a r)-i Sp '° Sq '° ! r = even) 
/ KP'9> r J- \ S p , S qtl (r = odd). 

In this way, we can implement the desired rotation Sz{9) up to Pauli byproducts X P Z 9 . Note that in the second 
term 

h(p,q,r)W(Z r \R) r+1 ) 

the desired rotation Sz(0) is not implemented (the trivial Z r is implemented, instead). However, this term can be 
treated as the usual error which can be corrected by the usual FT circuits. 

B. Simulation on the AKLT state with error 

Let us study what happens if an error occurs. We will see that not all physical errors on a physical qutrit can be 
linear CPTP errors in the correlation space. 

If a CPTP error occurs on the first physical qutrit, Eq. ([5]) becomes 

E E w (Qr(si, ■-, s r )...Q 2 {s ll s 2 )E^ Sl \R)) r +i ® -, s r ),g(s x , s r )) 



y-if N (\ L ),\R)) sj s 



' EE( E W(Q r ( Sl ,..., Sr )...Q 2 (s uS2 )E hSl \R)) r+1 



Vl /H " 3 P,1 (s 1 ,...,s r )eS pq 



+h(p, q, r)W{Z T -'E. h2 \R)) r +x) ® v(j>, «)■ 
This means that for fixed p and q, the map 

\R)(R\ -> E ^Qr(«i,...,ar)...Q2(ai,*a)£i,. 1 |ii)(iZ|4.iQ2(ai > *2)...Q3:(ai J ...,*r) 

( 8l ,.., 8r )es; i5 j 

+h(p, q, r) E Z r - 1 E J . 2 \R){R\E] 2 Z r - 1 
j 

is implemented in the correlation space. Note that 

E E El Si Ql( Sl ,s 2 )...Qt( Sl , s r )Q r (s u s r )...Q 2 (s u s 2 )E j , Sl + h(p, q, r) £ eX^-^ZT^E^ 

(si Oes*, j j 

E E 4.i + / " ; /'- ^ r ) E 4,2^,2 

(si,...,s r )eS; g j j 

= E E i T ™ i • E L E ^ + "■ r ) E E U E ^- 
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Here. 



I£* = {(si, s r ) S {0, 1, 2} xr \ (2, 2) | si = i and f(a u s r ) = p and g(si, s r ) 
For example, let us consider the case p = 1 and q = 0. As is shown in Appendix [El 



IT, 



1,0 I 



IT, 



1,0 1 



Therefore, 



V ^ lT r, sl pt P / KoV + Ej 4o%o (r = odd) 



For example, if w = 1 and 



^1 - ^M 9 ,^ ( 



|0) + |l)/n> |0> — 1 1> 



s/2 



V2 



1| + |2>(2| 



Then, 



^1,0 = y|o)(i| 
Si.i = y||i)<o| 



Wif*i F t F _ j |3lbV + ||i>UI (r = odd) 

2^2^ l J i,o I ' - i | T r.0|j + || \/ | + J/ f r = even), 10 



Sl j = l 



1,0 K ^ 3l"/\"l ^ 3 



+ il (r = even), 



which means that the map implemented in the correlation space is not linear CPTP. 



IV. INTUITIVE EXPLANATIONS 

So far, we have seen that the results of Refs. [TH, [lj| for the cluster state cannot be directly applied to other resource 
states, such as the one-dimensional AKLT state. Why the cluster state is so special? And why direct applications 
of the result for the cluster state to other resource states do not work? Although the complete answer to these 
questions is beyond the scope of the present paper, since the study of QCTN itself has not been fully developed (for 
example, no one knows the necessary and sufficient condition for tensor-network states to be universal resource states 
for measurement-based quantum computation), let us try to give some intuitive explanations here. 

Figure Q] illustrates the reason why all physical errors become CPTP maps in the correlation space of the one- 
dimensional cluster state. Let us first consider the ideal case (a) where there is no error. The actual protocol (a-1) is 
mathematically equivalent to the "input-output" picture (a-2) where the physical input state \ip) is teleported into the 
left-edge of the short- length chain (indicated in yellow) and finally the physical output state l^') is extracted from the 
right-edge of the short-length chain (indicated in yellow). Since both and \ip') are physical states, what is going 
on in the correlation space which maps the input state to the output state \if>') can be described by a linear CPTP 
operation. In other words, if we can describe measurement-based quantum computation with this "input-output" 
picture [l!|, the map implemented in the correlation space is guaranteed to be a linear CPTP operation (l9j . 

For the cluster state, this "input-output" picture also holds even if there is an error 19]: In the imperfect case, 
Fig. [1] (b), let us assume that the input state is degraded by an error and becomes a mixed state p. However, we can 
still consider a similar "input-output" picture (b-2), which corresponds to the actual protocol (b-1), and again the 
physical state p is mapped into another physical state p' , which means that what is going on in the correlation space 
which maps p to p 1 can be described by a linear CPTP operation. 
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Note that two special properties of the cluster state enable such an "input-output" picture. First, the one- 
dimensional cluster state can be decomposed into small pieces of one-dimensional cluster states by applying nearest- 
neighbour two-body unitary operations (i.e., CZ gates). As is shown in Fig. [T] (a-2) and (b-2), this property is 
necessary for allowing the "input-output" picture. Second, the number of qubits that are measured in order to imple- 
ment a specific gate does not depend on the measurement results. In other words, for the cluster state, a specific gate 
can be implemented up to Pauli byproducts at a fixed site irrespective of measurement results. Such a deterministic 
implementation at a fixed site is necessary for the deterministic (i.e., trace-preserving) "output" in the "input-output" 
picture, since in the method of Refs. [19(, the ensemble (mixture) of all measurement results are considered: If the 
site where the desired gate operation is completed depends on the measurement results, we cannot "extract" the same 
output state at a fixed site irrespective of measurement results as is shown in Fig. Q] (a-2) and (b-2). 

On the other hand, such an "input-output" picture seems to be impossible for the one-dimensional AKLT state, 
because of the following two reasons: First, as is shown in Fig.[5]left, no nearest-neighbour two-body unitary operation 
can decompose the one-dimensional AKLT chain into two chains due to the existence of the non- vanishing two-point 
correlation in the AKLT state. (If the one-dimensional AKLT chain can be decomposed into two chains by such a 
unitary, it contradicts to the well-known fact that the two-point correlation is non-vanishing in the AKLT state.) 
Second, we cannot deterministically implement a specific gate at a fixed site of the AKLT chain irrespective of 
measurement results [U, [H, |2{J. In short, the "input-output" picture seems to be impossible for the AKLT state. If 
we can no longer use the "input-output" picture, it is not unreasonable that we have some anomalous maps in the 
correlation space since the correlation space is not a physical space but an abstract mathematical space. 



(a) Ideal case 
(1) actual protocol 



OO 

input 



ooo 

output 

measurement 
(including feedforward) 

IV-') 

cx> 



(2) "input-output" picture 

I^L C + Z 



cz 
+ 



OO OOO 



— input output 

measurement 
mathematically (i nc | udin g feedforward) 
equivalent ^ . . 

W) C ^ 



OO 



(b) Imperfect case 
(1) actual protocol 

bo ooo- 

input output 
imperfect measurement 
(including feedforward) t 




imperfection 



OO 



(2) "input-output" picture 

"do- oo ? 0' 

input output 
— imperfect measurement 

(including feedforward) 
mathematically s. n' ^ 

OO 



equivalent 



FIG. 1: The "input-output" picture for the one-dimensional cluster state. 



V. ANOTHER WAY OF SIMULATING QUANTUM CIRCUITS 

As we have seen in the previous section, it is not always possible for general resource states to implement a specific 
gate at a fixed site of the resource state irrespective of measurement results. This fact prohibits a general resource 
state from allowing the "input-output" picture. One might think that if we abandon such a deterministic "output" 
at a fixed site for all measurement results, and if we just consider a specific history (trajectory) of measurement 
results, we might be able to avoid the emergence of non-CPTP errors. (Physically, this means that we project the 
system onto a pure state at every measurement step.) If the system is assumed to be error- free, this "trajectory 
method" is another standard way of simulating quantum circuits in the correlation space [3-0, E3 Hll . (Note that 
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non-zero two-point correlation between a and b 




FIG. 2: Left: The one-dimensional AKLT state cannot be decomposed into short chains. Right: The tensor network for the 
resource state in Ref. [3l| . 



in this trajectory method, correct unitary operators can be implemented in the correlation space if there is no error, 
although what we are physically doing are projections, i.e., non-trace-preserving operations). 

However, we here show that such a natural another way of simulating quantum circuits does neither work if d > 3. 
In other words, we can show the following theorem. 



A. Theorem 



Theorem: If d > 3, there exists a single-qudit CPTP error £ which has the following property: assume that £ is 
applied on a single physical qudit of \*f?(L, R))^ . If the measurement Aie.<p, Eq. ([2]), is performed on that affected 
qudit, a non-TP operation is implemented in the correlation space. 

Proof: In order to show Theorem, let us assume that 

There is no such £. (8) 

We will see that this assumption leads to the contradiction that d < 2. 
First, let us consider the state 

(7® JV - 1 ®C/ 1 „ 2 )|*(L,i2))f, (9) 

where 



U a ^ b = \a) (b\ + \b) (a\+I- \a) (a\ - \b) (b\ 

is the unitary error which exchanges \a) and \b), and I is the identity operator on a single qudit. In Eq. (J9j) , the error 
is applied on the first physical qudit of \^f(L, R))i ■ If we do the measurement M.e,<j> on the first physical qudit 
of Eq. ©, and if the measurement result is |2), Eq. © becomes 

-—-^==—\\A[l]\\ J2 .» E MA[k N }..A[k 2 ] 1 ^\R)\k N ,...,h)®\2). (10) 
VfN-i(\L),A[l]\R)) k % k ^ Q \\A[1]\\ 

In other words, the operation 

|R> "* PM |fi> 

is implemented in the correlation space. By the assumption Eq. (|8|), this operation should work as a TP operation in 
the correlation space. Therefore, 

At W A W i. (ii) 



\\A[1] 

By taking -q = ||A[1]|| 2 , 

A^[l]A[l]=rjI. (12) 



Second, let us consider the measurement A4e,<f> on the first physical qudit of 



JV 

1 ) 



where s £ {0, 1, d — 1}, 



V ^J2e-^\p)(p\ 

p=0 

is a unitary phase error, and w = 2it / d. If the measurement result is \otg 

e -2fe« cos °_ A [2] + e -'(0+^) sin ^[i]) /^y 
is implemented in the correlation space, where 

yfy= e~ 2isu cos -A[2] + e'^ 8 ^ sin -A[l] . 
By the assumption Eq. ([5]). this should work as a TP operation in the correlation space. Therefore, 

1 I = cos 2 -At[2]A[2] + sin 2 -^t[l]A[l] + - smO^-^-*^ A^[2]A[1] + e^^A* [l]A[2]j . 
By the assumption Eq. ([3]), 

Al[2]A[2]=£I, 

where £ = ||^4[2]|| 2 . Furthermore, as we have shown, = r/I (Eq. (I12[) V Therefore, 

yj = e -*W-«-)A+[2]A[l] + e^^Ut^Ap], 

where 

2 



7 



sin i 



cos 2 77 sin 2 — 

2 ' 2 



Finally, let us consider the measurement Me,4> on the first physical qudit of 
where t £ {0, 1, d — 1}. If the measurement result is \ag^), 



V**" C o S -A[l] + e -^- 2it « S in -A[2]j /V* 



is implemented in the correlation space, where 



V~8 = 



»-**" cos -All] + e -^- 2jt " sin -A{2] 
2 1 1 2 1 1 



By the assumption Eq. ©, this should also work as a TP operation in the correlation space. Therefore, 

g'j _ e ^+M^t[ 2 ]A[l] + e - i ^+' a 'Ut[l]A[2], 

where 



5' 



sm( 



r i- ■ 2 u 2 U 

o — t sin 77 cos - 

s 2 ' 2 



From Eqs. {13]) and flUJ), 



e -2i(<j>-su>) _ e 2i(<j>+tuj) 



tf[2]A[l], 
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where 

e = e _,; ^ _sw - ) 7' — e i (<H-* a; ) ( 5 / , 

Let us assume that 

e -2i(4>-suj) e 2i(<p+tu>) i g 

Then, 

e'I = A*[2]A[1], 

where 



g— 2i(0— sw) g2i(0-fio;) 

If e' = 0, A^[2]A[1] = 0, which means A[l] = since A[2] is unitary up to a constant (assumption Eq. Therefore, 
e' 7^ 0. In this case, A[l] = e"A[2] for certain e" ^ 0, since A [2] is unitary up to a constant [3J]. Hence 

e -2j(0-sw) _ e 2i(0+tw) _ g 

This means 

2cf) + (t - s)uj = r 8>t TT, (15) 
where r 8 ,t G {0, 1,2,3, ...}. Let us take t = s = 0. Then, Eq. (fTS"|) gives 

^ = ro,o| (r , e {0,1, 2, ...}). 

Let us take s = 1, t = 0. Then, Eq. (fT5j) gives 

^=5 +ri '°f (n,o S {0,1,2,...}). 

In order to satisfy these two equations at the same time, there must exist ro,o and r^o such that 

7T 7T IT 

If ro,o = J"i,o> then = 1/d which means d = oo. Therefore ro,o 7^ ^l.o- Then we have 

2 

d = < 2, 

^o,o - n,o 

which is the contradiction. ■ 

One might think that if we rewrite the post-measurement state Eq. (|10p as 

V/»-,(|iU]W) IWWI1 £ -E^M-^Np^iWI*-. ..,*> « P) 

and redefine the operation implemented in the correlation space as 

|i?) ^RllM |i?> ' 

the TP-ness is recovered in the correlation space. However, in this case, the non-lineally appears unless 

A f [l],4[l] oc I, 

and therefore if we require the linearity in the correlation space, we obtain the same contradiction. 

In short, if d > 3 not all physical errors on physical qudits appear as linear CPTP errors in the correlation space of 



pure matrix product states [35 1. 
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B. Examples 



Let us consider some concrete examples which give intuitive understandings of the above theorem. 

In Ref. [29j, it was shown that universal single-qubit unitary rotation is possible in the correlation space of the 
one-dimensional AKLT chain, which is a ground state of a gapped two-body nearest-neighbour spin-1 Hamiltonian 
(hence d = 3). The matrix product representation of the one-dimensional AKLT chain is given by 



for certain basis {|0), |1), |2)} [H 



A[0] 
A[l] 
A[2] 

The measurement 



Me 



,tt/2 = i COS -| 



i sin — 11) 



xz, 
z 



,sm-|0) 



l cos 



|lD,»} 



on a single physical qutrit implements Xe~ lZ9 / 2 , XZe~ %Z6 1 2 , or Z, respectively. According to Theorem in the previous 
section, not all physical errors can be linear CPTP maps in the correlation space since d > 3. In fact, let us consider 
the single qutrit unitary error 



U = \2)- 



(II 



10} 



Hi 



V2 



V2 



<2| 



|Q)-|i) (Q|-(i| 

V2 V2 



If the measurement is performed after the error U and if we obtain the result |2), the operation |1)(0| is 

implemented in the correlation space. Obviously, it is not TP. The same result is obtained for the slightly modified 
version of the one-dimensional AKLT chain jij, where 

A[0] = X, 
A[l] = XZ, 
A[2] = H. 

In Ref. [31], it was shown that the unique ground state of a gapped two-body nearest- neighbour spin-3/2 Hamiltonian 
with the AKLT and exchange interactions on the two-dimensional octagonal lattice is a universal resource state for 
the measurement- based quantum computation. The state is defined by the following tensor network |3lj : 

= |1><0|®<1|, 



At 


r 3 i 

[ + 2\ 


At 


■ 3 
.~ 2. 


At 


l + 2\ 


At 


■ 1- 

."" 2. 


D 




D 


- 3 
.~ 2. 


D 


r ii 

[ + 2\ 


D 


■ 1- 
.~ 2. 



(1| + |1)<0|®<0| 
(z®(O|-|O)(l|0(l 



|0)(1|®<0|, 

— \={z® 

V3V 

1 

|0)(1|, 

-|i)(o|, 



(see Fig. fright). ^4^'s are defined in the same way. Each horizontal line works as a single-qubit wire. Two nearest- 
neighbour horizontal chains are decoupled by measuring sites B in the z-basis. Before starting the computation, the 
filtering operation {F, F' }, where 



F 
F' 



1 

- 



|3/2) <3/2|- 
|3/2><3/2| 



3/2)<-3/2| + ]l/2>(l/2| 
- 3/2><-3/2|, 



l/2><-l/2| 
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is applied on each site A. Let us assume that the filtering is succeeded (i.e., F is realized) and a site B is projected 
onto 1 3/2). Then, the measurement 

= {^=(11/2) ± e*| - 3/2)), | - 1/2), | + 3/2)} 

implements ZXe lZ4>/2 , le tZ * /2 , or Z, respectively. (The result |3/2) does not occur.) If the error which exchanges 
|l/2) and | — 1/2) occurs and if the result of the measurement M.^/ 2 ^ is | — 1/2), the operation |1)(0|, which is not 
TP, is implemented in the correlation space. 



VI. CONCLUSION 



In this paper, we have studied how physical errors on a physical qudit appear in the correlation space of general 
resource states. We have shown that the results [3, [l9[ for the cluster state cannot be directly applied to general 
resource states, such the AKLT state. We have also shown that if d > 3 not all physical errors can be linear CPTP 
errors in the correlation space of pure matrix product states. These results suggest that the application of the theories 
of fault-tolerant quantum circuits to the correlation space of general resource states is not so straightforward. 
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Appendix A: Tricluster state without error 

Let us consider the tricluster state [23], where d = 6 and 



A[0] = 


]+)<o 


A[l] = 


|->(1 


A[2] = 


|->(o 


A[3] = 


]+)<! 


A[4] = 


]+)<! 


A[5] = 


|-)<o 



The measurement 



o) = 




JO) 


+ e u 


Jl>) 


l) = 




JO) 


-e lt 


Jl>) 


2) = 




J2) 


+ e u 


J3)) 


3) = 


U 


J2) 


-e u 


'|3>) 


4) = 




J4) 


+ e~ 


»|5>) 


5) = 




J4) 


— e~ 


»|5>) 
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is performed in order to implement 



\0o) 


J(0) 




XJ(0) 


l&) 


ZJ{6) 


l*3> 


ZXJ{6) 


1^4) 


ZJ{6) 


l<?5> 


ZXJ{6) 



First, we measure the first physical qudit in the {|0 S1 )} basis. Then we obtain 

5 



1 



2f N (\L),\R)) 



W(X^Z^\J(6)\R)) 2 ® \6 S1 ){6 S1 \® m(ai), 



si=0 



where 



p(s) = 
Q(s) = 

If we trace out the first physical qudit, we obtain 



S s ,i + <5 S ,3 + S 8l 5 

S s .2 + S s .3 + SsA + S a R. 



Second, if we measure the second physical qudit in the basis 

l-i, A.) = -^(|0> + (-l)^)e**(- 1 ) P( - 



\si,<pi) = 
\si,(h) = 

\si,(f>4) = 
\S1A5) = 

and trace out the measured qudit, we obtain 

5 5 



V2' 

_2_Q > - (_i)«(«x) e ^<- 1 >" <si 

-^=(|2) + (-l)^) e ^- 1 ) P<ai 

-i=(|2)-(-l)»(' 1 )e i *(- 1 ) ,,( * 1 
v2 ^ 

_L(| 4) + ( _l )9 ( S i) e #(-i) p(si 
a/2 ^ 



(Al) 



1 



2*f N (\L),\R)) 



2*f N (\L),\R)) 



W(^ p(s2) ^ 9(;S2) ^ < ' (;Sl) ^((-l) p(;Sl) 0)^ p(si) ^ (;Sl) ^(6 l )|i?))3 ® m(«i) ® m(s 2 ) 



S 2 =0 Sl=0 

5 5 



m(s 2 )- 



S 2 =0 Sl=0 
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Third, if we measure the third physical qudit in the basis 



|si,s 2 ,?7o) 
|si,S2,»7x) 

|S1,S2,»73) 
|S1,S2,»?4) 
|S1,S2,»?5) 

\si,s 2 ,rj 6 ) 



V2 
1 

71 
i 

71 
i 

71 

U 
U 



|o> 

|0) 
|2> 
|2) 
|4> 

|4) 



. 1 \p(Bi)+q(s 2 ) e ir,(-iy<-^ 



.l)p(si)+q(s 2 ) e iri(-l) Pls2 



.iyM+q{s 2 ) e ir,(-ir<-^ 



_l)P(si)+9(s2) e ir,(-l)f< s 2 



and trace out the measured third qudit, we obtain 

5 5 5 

EEE 



i 



(A2) 



2 3 /Ar(|L),|i?)) 

W ( XP (s 3 ) Zq ( S3 ) X p( Sl )+ q (s 2 ) J^ l Y(s2)^ X p(s2) Z q{s2) Z p{s 1 ) J(0)\R)) 4 (g) m ( Sl ) O m ( S2 ) <g) m ( S3 ) 

5 5 5 



1 



23f N (\L),\R}) 
1 



EEE W ( XP[S3) zq{S3) ZP(S2) J(vW)JW)\R)h 

® m(si) ® m(s2) <8> m(s3) 



S3— S2— Si— 

5 5 



23 /Ar (|L),|i?))^ 0s f ro 



J2 E W / (^ p(S3) ^ (S3) ^ p(S2) ^('7)^(0)>/(^)| J R))4®m( S 2)®m( S 3). 



Appendix B: Tricluster state with error 



Let us assume that a CPTP error occurs on the first physical qudit. If we measure the first physical qudit in the 
{|0 S i)} basis, we obtain 

1 5 

By tracing out the first physical qudit, we obtain 

Second, if we measure the second physical qudit in the basis Eq. (|A1|) and trace out the measured qudit, we obtain 

5 5 

2f (lL) lR)) E T,zZ w ( XP(S2)zq{S2)xq{sl)j ^- 1 ^ 

jn{\ ), I /J S2=0 si=0 ^ 



Third, if we measure the third physical qudit in the basis Eq. (| A2|) and trace out the measured third qudit, we 
obtain 



1 



5 5 5 



E E EE 



2 2 f N (\L)AR)) ^ 

Vt/(x p ( S3 )Z9( S3 )XP( sl ) + «( S2 ) j((-i)p(«) r/ )x p ( S2 )Z' ? ( S2 )X <? ( S1 ) J((-l) p ( sl V)^j,s 1 |-R»4 ® m(si) ® m(s 2 ) ® m(s 3 ) 

5 5 5 

2 2 f N (\L),\R)) E E E E 

jjv vi / i n S3=0 S2=0 Si=0 J 
jy^Jz*)^)^'!) J[ri)X q ^ J((-1) P(S1 tyj-E*,., |.R))4 ® m(si) ® m(s 2 ) <8 m(s 3 ). 



If we trace out the first record, 

1 5 5 5 

22f m m E E ^ E ^P' 83 '^*^* 2 '^^^!^'''^!!-!)*''^^,..^})! ® m( S2 ) « 

JN{\ },\ )) S3=0s2=0si=0 j 

Thus the map 

5 

\R)(R\ -> E E^ ^)^ 1 ^^ -1 )^ 1 ^)^!^^!^! ^((-l)^' ^^" ^^)^' , 

«i=0 j 

which is obviously TP, is implemented. 



Appendix C: Calculation of \Up 

Let us define 



U, 



p, q = |(si,-,Sr) e {0, l,2} xr /(si, ...,s r ) = p and g(s 1 , ...,s r ) = g|. 



First, 



l^ol 


= 3 


|c/ 2 il 


= 2 


|C/i 2 .ol 


= 2 




= 2. 



Second, 

l^p.gl = l^pffil,(jl + l^p,g©ll + l^pei,g©il 

for all r. Therefore, 

l^il-l^ol = l^l + l^l + l^oVl^oVl^T 1 !-!^,! 1 
= l^o'l-I^I 1 ! 
= -(Kl'l-Ko 1 !) 
= (-ir 2 (l^,il-l^i 2 ol) 

= 



and 



which mean 



for all r. 
Note that 



Kil-Kil = l^l + l^l + l^o 1 !-!^! 1 !-!^ 1 !-!^ 1 ! 

= -(I^T 1 !-!^,! 1 !) 

= (-ir 2 (l^o 2 ,il-l^i 2 ,il) 
= o, 



It^oM^il = 3|t/ 7 i |-(2|C/o r ,T i | + l^,ol) 

= l^l'l-l^o 1 ! 

= -d^oVl^T 1 !) 

= (-ir 2 (i^o 2 oi-i^o 2 ii) 
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and 



Therefore, 



Kol+3|£/ r ,il 



= 3|[/ 7 1 |+3(2|C/ ^ 1 | + |C/ 7 1 |) 

= si^i + qi^y 1 ! 

= y- 2 (\ul a \ + '3\uL\) 

= 3 r . 



= -(3 r -(-iy) 
\u£, \ = l(3 r + 3(-iy). 



Appendix D: Calculation of \Sp } , 



Let us define 



S r piq = {(«!, s r ) G {0, 1, 2} xr \ (2, 2) I f( Sl , s r ) = p and s r ) = q} 



First, 



Second, 



Therefore, 



and 



Hence 



...,2) 


f(si 


^0,0 


= 2 


|So 2 il 


= 2 


|S?,ol 


= 2 


\sh\ 


= 2. 



I^p.gl — l^p©l,J + l^p©l,g©ll + l^p,g©ll- 

l^,o|-|^o,il = l^l + l^l + l^iVl^l-l^oVl^o 1 ! 

= -d^o 1 !-!^,! 1 !) 
= (-i) r - 2 (|sg,ol-|sg,il) 

= 0. 

l^,ol-l^,il = l^l + l^l + l^l-I^I 1 !-!^ 1 !-!^ 1 ! 
= -d^o 1 !-!^,! 1 !) 

= (-ir'd^.oi-i^i) 

= 0. 



\Sq,o\ 



3 r-l _ (_iy 



r-li 



+ |5 r ,o 



s 7 - 1 - (-1) 



r-l 



3 2 -(-l) 2 



2 2 
3'- 1 + ... + 3 2 (-ir- 1 + ... + (-l) 2 
2 2 

3 r -3 2 (-l) r -(-l) 2 



K \ 
+ 2 



+ 



+ 2. 



and 



\s[, \ 



y- 1 + (-1)' 



r-l| 



+ \s[; 



r-l 



3 2 + (-l) : 



2 - ' 2 ' |5l >° 

S^i + .-. + S 2 , (-l)''- 1 + ... + (-l) 2 



2 + 2 

3 r -3 2 (-ir-(-i) 2 , 2 



+ 2 



It is easy to see that 



Appendix E: Calculation of \Tpf q \ 



IT''' 1 ! — ITT''" 1 I 

\ ± p,q\ ~ \ u p®l,q®l\ 

|7"T,2 1 I or— 1 I 

\ ± p,q\ ~ l°p,g©lh 



where 



Since 



U P,q = {( s i'-> s r) G {0, l,2} xr /(si,...,s r ) =p and #(si, s r ) 



1^1 = l^il = 7( 3r_1 -(- 1 ) r " 1 ) 



1^1 



r-l 



we obtain |T r ;°| = \T^\ and 



i(3 r - 1 -3 2 + (-l) I '- 1 -(-l) 2 ) + 2, 



-1 3 2 1 „ 
1 1 2 

2 4 4 

r-l 



l-(-l) 

2 

r (r = odd) 



= U 



(r = even). 



Since 



\T$\ 



I^.T 1 ! 



i (3 f - 1 + 3(-ir 1 ) 
l 

4 

3 r-l _ 32 



(3*- 1 - (-if- 1 ) 

_ X )r-l _ (_ 1)2 



+ 2, 



we obtain 



l^'oVKol = (-!)' 



i^i-i^'oi = ^ + (-ir 1 -J-2 
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and 



Since 



we obtain 



and 



= 0. 



\T r J\ = \U^ 1 1 \ = \(3 r - 1 -(-l) r - 1 ) 

= i^,o 1 i = ^ 3r " 1 + 3 (- 1 ) r " 1 ) 



3 r-l_ 3 2 (-l)r-l _ (_!)2 



1^1 -1^1= -(-If" 1 , 



Tffl - iTTfl = - ---2 



3 2 1 



-i,i I l- 1 1,1 



4 4 
0, 



i^i-i^i - T + Hr-7-2 



^ 2 i 

4 V ; 4 
(-If- 1 . 
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